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We present and discuss the results of calculations of surface relaxations and rumplings for the 
(001) and (011) surfaces of BaTiOs and PbTiOa, using a hybrid B3PW description of exchange and 
correlation. On the (001) surfaces, we consider both AO (A = Ba or Pb) and Ti02 terminations. 
In the former case, the surface AO layer is found to relax inward for both materials, while outward 
relaxations of all atoms in the second layer are found at both kinds of (001) terminations and 
for both materials. The surface relaxation energies of BaO and Ti02 terminations on BaTiOa 
(001) are found to be comparable, as are those of PbO and Ti02 on PbTiOa (001), although 
in both cases the relaxation energy is slightly larger for the Ti02 termination. As for the (011) 
surfaces, we consider three types of surfaces, terminating on a TiO layer, a Ba or Pb layer, or an 
O layer. Here, the relaxation energies are much larger for the TiO-terminated than for the Ba or 
Pb-terminated surfaces. The relaxed surface energy for the O-terminated surface is about the same 
as the corresponding average of the TiO and Pb-terminated surfaces on PbTi03 , but much less than 
the average of the TiO and Ba-terminated surfaces on BaTiOs. We predict a considerable increase 
of the Ti-O chemical bond covalency near the BaTiOs and PbTiOs (011) surface as compared to 
both the bulk and the (001) surface. 

PACS numbers: 68.35.bt, 68.35.Md, 68.47.Gh 



I. INTRODUCTION 

Thin films of ABO3 perovskite ferroelectrics play an 
important role in numerous microelectronic, catalytic, 
and other high-technology applications, and are fre- 
quently used as substrates for growth of other materi- 
als such as cuprate superconductors*^ Therefore, it is 
not surprising that a large number of ab initio quan- 
tum mechanical calculations ^^^^ i 10 ' 11 ' 12 ' 13 as well 
as several classical shell-model (SM) studies ; 14 ' 15 ! 16 have 
dealt with the atomic and electronic structure of the 
(001) surface of BaTi0 3 , PbTi0 3 , and SrTi0 3 crystals. 
In order to study the dependence of the surface relax- 
ation properties on the exchange-correlation functionals 
and the type of basis (localized vs. plane-wave) used in 
the calculations, a detailed comparative study of SrTi0 3 
(001) surfaces based on ten different quantum-mechanical 
tcchnique a 17 i 18 was recently performed. 

Due to intensive development and progressive minia- 
turization of electronic devices, the surface structure as 
well as the electronic properties of the AB0 3 perovskite 
thin films have been extensively studied experimentally 
during the last years. The SrTi0 3 (001) surface struc- 
ture has been analyzed by means of low-energy elec- 
tron diffraction (LEED), 19 reflection high-energy elec- 
tron diffraction (RHEED), 2 — X-ray photoelectron spec- 
troscopy (XPS), ultraviolet electron spectroscopy (UPS), 
medium-energy ion scattering (MEIS)f 21 - and surface X- 
ray diffraction (SXRD). 22 Nevertheless it is important to 
note that the LEED 1 ^ and RHEED 20 experiments con- 
tradict each other in the sign (contraction or expansion) 
of the interplanar distance between top metal atom and 
the second crystal layer for the SrO-terminated SrTi0 3 
(001) surface. The most recent experimental studies 



on the SrTi0 3 surfaces include a combination of XPS, 
LEED, and time-of- flight scattering and recoil spectrom- 
etry (TOF-SARS), 23 as well as metastable impact elec- 
tron spectroscopy (MIES).— In these recent studies, well- 
resolved lxl LEED patterns were obtained for the TiOV 
terminated SrTi0 3 (001) surface. Simulations of the 
TOF-SARS azimuthal scans indicate that the O atoms 
are situated 0.1 A above the Ti layer (surface plane) in 
the case of the Ti02-terminated SrTi0 3 (001) surface. 

While the (001) surfaces of SrTi0 3 , BaTi0 3 and 
PbTi0 3 have been extensively studied, much less is 
known about the (011) surfaces. The scarcity of informa- 
tion about these surfaces is likely due to the polar charac- 
ter of the (011) orientation. (011) terminations of SrTi0 3 
have frequently been observed, but efforts towards the 
precise characterization of their atomic-scale structure 
and corresponding electronic properties has only be- 
gun in the last decade, specifically using atomic-force 
microscopy,^ 5 " scanning tunneling microscopy (STM), 
Auger spectroscopy, and low-energy electron-diffraction 
(LEED)^ methods. 

To the best of our knowledge, very few ab initio stud- 
ies of perovskite (011) surfaces exist. The first ab initio 
study of the electronic and atomic structures of several 
(lxl) terminations of the (011) polar orientation of the 
SrTi0 3 surface was performed by Bottin et al£- One 
year later Heifets et al£^ performed very comprehensive 
ab initio Hartree-Fock calculations for four possible ter- 
minations (TiO, Sr, and two kinds of O terminations) 
of the SrTi0 3 (011) surface. Recently Heifets et al. 29 
performed ab initio density-functional calculations of the 
atomic structure and charge redistribution for different 
terminations of the BaZr0 3 (011) surfaces. However, 
despite the high technological potential of BaTi0 3 and 
PbTi0 3 , we are unaware of any previous ab initio cal- 
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culations performed for the BaTiOs and PbTi03 (Oil) 
surfaces. In this study, therefore, we have investigated 
the (Oil) as well as the (001) surfaces of BaTiOs and 
PbTiC>3, with an emphasis on the effect of the surface 
relaxation and rumpling, surface energies, and the charge 
redistributions and changes in bond strength that occur 
at the surface. 



II. PRELIMINARIES 

A. Computational method 

We carry out first-principles calculations in the frame- 
work of density-functional theory (DFT) using the 
CRYSTAL computer code^ Unlike the plane- wave codes 
employed in many previous studieS ) 31 i 32 CRYSTAL uses 
localized Gaussian-type basis sets. In our calculations, 
we adopted the basis sets developed for BaTi03 and 
PbTi03 in Ref. [HI]. Our calculations were performed 
using the hybrid exchange-correlation B3PW functional 
involving a mixture of non-local Fock exact exchange, 
LDA exchange, and Becke's gradient corrected exchange 
functional? 3 ^ combined with the non-local gradient cor- 
rected correlation potential of Perdew and Wang i 35 i 36 i 37 
We chose the hybrid B3PW functional for our current 
study because it yields excellent results for the SrTi03, 
BaTi03, and PbTiOs bulk lattice constant and bulk 
modulus j 17 ' 33 

The reciprocal-space integration was performed by 
sampling the Brillouin zone with an 8x8x8 Pack- 
Monkhorst meshf 3 ^ which provides a balanced summa- 
tion in direct and reciprocal spaces. To achieve high 
accuracy, large enough tolerances of 7, 8, 7, 7, and 
14 were chosen for the dimensionless Coulomb overlap, 
Coulomb penetration, exchange overlap, first exchange 
pseudo-overlap, and second exchange pseudo-overlap pa- 
rameters, respectively^ 

An advantage of the CRYSTAL code is that it treats 
isolated 2D slabs, without any artificial periodicity in the 
z direction perpendicular to the surface, as commonly 
employed in most previous surface band-structure cal- 
culations (e.g., Ref. 0). In the present ab initio inves- 
tigation, we have studied several isolated periodic two- 
dimensional slabs of cubic BaTi03 and PbTiOs crystals 
containing 7 planes of atoms. 

B. Surface geometries 

The BaTiC-3 and PbTi0 3 (001) surfaces were modeled 
using symmetric (with respect to the mirror plane) slabs 
consisting of seven alternating Ti02 and BaO or PbO 
layers, respectively. One of these slabs was terminated 
by BaO planes for the BaTiOs crystal (or PbO planes 
for PbTiOs) and consisted of a supercell containing 17 
atoms. The second slab was terminated by Ti02 planes 
for both materials and consisted of a supercell containing 




FIG. 1: (Color online.) Side view of a Ti02-terminatcd 
BaTiOs (001) surface with the definitions of the surface rum- 
pling s and the near-surface interplanar distances di2 and t?23 , 
respectively. 



18 atoms. These slabs are non-stoichiometric, with unit 
cell formulae Ba4Ti30io or Pb4Ti30io, and Ba3Ti40n 
or Pb 3 Ti 4 0n for BaTi0 3 and PbTi0 3 perovskites, re- 
spectively. These two (BaO or PbO and Ti02) termi- 
nations are the only two possible flat and dense (001) 
surfaces for the BaTiOs or PbTiOs perovskite lattice 
structure. The sequence of layers at the Ti02-terminated 
(001) surface of BaTi03 is illustrated in Fig. 1. 

Unlike the (001) cleavage of BaTi0 3 or PbTi0 3 , which 
naturally gives rise to non-polar BaO (or PbO) and Ti02 
terminations, a naive cleavage of BaTiOs or PbTiOs to 
create (011) surfaces leads to the formation of polar sur- 
faces. For example, the stacking of the BaTiOs crystal 
along the (011) direction consists of alternating planes 
of O2 and BaTiO units having nominal charges of — 4e 
and +4e respectively, assuming 2 ~, Ti 4+ , and Ba 2+ 
constituents. (Henceforth we shall use BaTiOs for pre- 
sentation purposes, but everything that is said will apply 
equally to the PbTiOs case.) Thus, a simple cleavage 
leads to 02-terminated and BaTiO-terminated (011) sur- 
faces that are polar and have nominal surface charges of 
— 2 e and +2 e per surface cell respectively. These are 
shown as the top and bottom surfaces in Fig. 2(a) re- 
spectively. If uncompensated, the surface charge would 
lead to an infinite electrostatic cleavage energy. In re- 
ality, the polar surfaces would probably become metal- 
lic in order to remain neutral, but in view of the large 
electronic gaps in the perovskites, such metallic surfaces 
would presumably be unfavorable. Thus, we may expect 
rather generally that such polar crystal terminations are 
relatively unstable in this class of materials* 3 - 

On the other hand, if the cleavage occurs in such a way 
as to leave a half layer of O2 units on each surface, we ob- 
tain the non-polar surface structure shown in Fig. 2(b). 
Every other surface O atom has been removed, and the 
remaining O atoms occupy the same sites as in the bulk 
structure. We shall refer to this as the "O-terminated" 
(011) surface, in distinction to the "02-terminated" po- 
lar surface already discussed in Fig. 2(a). The non-polar 
nature of the O-terminated surface can be confirmed by 
noting that the 7-layer 15-atom BasTisOg slab shown in 
Fig. 2(b), which has two O-terminated surfaces, is neu- 
tral. It is also possible to make non-polar TiO-terminated 
and Ba-terminated surfaces, as shown in Figs. 2(c) and 
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FIG. 2: (Color online.) Side views of slab geometries used to 
study BaTiOs (Oil) surfaces, (a) Stoichiometric 8-layer slab 
with 02-terminated and BaTiO-terminated surfaces at top 
and bottom respectively, (b) 7-layer slab with O-terminated 
surfaces, (c) 7-layer slab with TiO-terminated surfaces, (d) 
7-layer slab with Ba-terminated surfaces. 

(d), respectively. This is accomplished by splitting a 
BaTiO layer during cleavage, instead of splitting an O2 
layer. For the TiO and Ba-terminated surfaces, we use 
7-layer slabs having composition Ba2Ti40io (16 atoms) 
and Ba4Ti 2 0s (14 atoms) as shown in Fig. 2(c-d), respec- 
tively. These are again neutral, showing that the surfaces 
are non-polar (even though they no longer have precisely 
the bulk BaTiC"3 stoichiometry) . 

III. RESULTS OF CALCULATIONS 
A. BaTi0 3 and PbTi0 3 (001) surface structure 

In the present calculations of the BaTiOs and PbTiC>3 
(001) surface atomic structure, we allowed the atoms lo- 
cated in the two outermost surface layers to relax along 
the z-axis (the forces along the x and y-axes are zero by 
symmetry). Here we use the term "layer" to refer to a 
BaO, PbO, or TiC>2 plane, so that there are two layers 
per stacked unit cell. For example, on the BaO or PbO- 
terminated surfaces, the top layer is BaO or PbO and 
the second layer is Ti02; displacements of the third-layer 
atoms were found to be negligibly small in our calcula- 
tions and thus were neglected. 

The calculated atomic displacements for the Ti02 and 
BaO-terminated (001) surfaces of BaTiOs, and for the 
Ti0 2 and PbO-terminated (001) surfaces of PbTi0 3 , are 
presented in Table I. For BaTiOs (001), comparisons 
are also provided with the surface atomic displacements 
obtained by Padilla and Vanderbilt 4 using plane-wave 
DFT methods in the local-density approximation (LDA) , 
and by Heifets et al. using a classical shell-model (SM) 



approach^ Similarly, for PbTi0 3 (001), Table I shows 
comparisons with the plane-wave LDA calculations of 
Meyer and Vanderbilt. 4 The relaxation of the surface 
metal atoms in both the BaTiOs and PbTi03 surfaces 
is much larger than that of the oxygen ions, leading to 
a considerable surface rumpling, which we quantify via a 
parameter s defined as the relative displacement of the 
oxygen with respect to the metal atom in a given layer. 
The surface rumpling and relative displacements of three 
near-surface planes are presented in Table II. According 
to our calculations, atoms of the first surface layer re- 
lax inwards (i.e., towards the bulk) for BaO and PbO 
terminations of both materials. Our calculations are in 
a qualitative agreement with the ab initio calculations 
performed by Padilla and Vanderbilt 4 for BaTiOs, and 
by Meyer and Vanderbilt^ for PbTi03. However, the 
predictions of the SM calculation disagree with the first- 
principles calculations; the SM predicts that the first- 
layer oxygen ions relax outward on the BaO-terminated 
BaTiOs (001) surface^ rather than inwards. However, 
the magnitudes of the displacements are relatively small 
(—0.63% of the lattice constant ao in this study and 
1.00% of ao in the SM calculations) which may be close 
to the error bar of the classical shell model. Outward re- 
laxations of all atoms in the second layer are found at 
both (001) terminations of the BaTi0 3 and PbTi0 3 sur- 
faces. From Table I, we can conclude that the magnitudes 
of the atomic displacements calculated using different ab 
initio methods and using the classical shell model are in 
a reasonable agreement. 

In order to compare the calculated surface structures 
further with experimental results, the surface rumpling s 
and the changes in interlayer distances Adi2 and Ag?23, 
as defined in Fig. 1, are presented in Table II. Our calcu- 
lations of the interlayer distances are based on the posi- 
tions of relaxed metal ions, which are known to be much 
stronger electron scatters that the oxygen ionsJ^ 

For BaTiOs (001), the rumpling of Ti02-terminated 
surface is predicted to exceed that of BaO-terminated 
surface by a factor of two. This finding is in line with 
the values of surface rumpling reported by Padilla and 
Vanderbilt, 4 In contrast, PbTi03 demonstrates practi- 
cally the same rumpling for both terminations. From 
Table II one can see that qualitative agreement between 
all theoretical methods is obtained. In particular, the re- 
laxed (001) surface structure shows a reduction of inter- 
layer distance Adi2 and an expansion of Ag?23 according 
to all ab initio and shell- model results. 

As for experimental confirmation of these results, 
we are unfortunately unaware of experimental measure- 
ments of Adi2 and Ac?23 for the BaTiOs and PbTi03 
(001) surfaces. Moreover, for the case of the SrO- 
terminated SrTi03 (001) surface, existing LEED 19 and 
RHEED 20 experiments actually contradict each other re- 
garding the sign of Adi2 . In view of the absence of clear 
experimental determinations of these parameters, there- 
fore, the first-principles calculations are a particularly 
important tool for understanding the surface properties. 



4 



TABLE I: Vertical atomic relaxations (in percent of bulk lattice constant) for BaTi03 and PbTiC>3 (001) surfaces. Positive sign 
corresponds to outward atomic displacement. 'SM' indicates shell-model calculation of Ref. [l6l |: 'LDA' are previous calculations 
of Ref. H and H for BaTi0 3 and PbTiQ 3 respectively. 
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B. BaTi0 3 and PbTi0 3 (Oil) surface structure 

To our knowledge, we have performed the first ab 
initio calculations for BaTiC>3 and PbTiC>3 (011) sur- 
faces. We have studied the Ti02-terminated, BaO or 
PbO-terminated, and O-terminated surfaces illustrated 
in Fig. 2(c), (d), and (b), respectively. The computed 
surface atomic relaxations are reported in Table III. 

Focusing first on the BaTiC>3 surfaces, we find that 
the Ti ions in the outermost layer of the TiO-terminated 
surface move inwards (towards the bulk) by 0.0786 ao, 
whereas the O ions in the outermost layer move outwards 
by a 0.0261 ao- The Ba atoms in the top layer of the 
Ba-terminated surface of Fig. 2(d) and the O atoms in 
the outermost layer of O-terminated surface of Fig. 2(b) 
move inwards by 0.0867 a and 0.0540a , respectively. 
The agreement between our ab initio B3PW and the clas- 
sical SM calculations is satisfactory for all three of these 
surface terminations. In particular, the directions of the 
displacements of first and second-layer atoms coincide for 
all three terminations. This indicates that classical SM 
calculations with a proper parameterization can serve as 
a useful initial approximation for modeling the atomic 



TABLE II: Calculated surface rumpling s and interlayer dis- 
placements Adij (in percent of bulk lattice constant) for near- 
surface planes on the BaO or PbO and TiC>2-terminated (001) 
surfaces of BaTiQ 3 and PbTi0 3 . 
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structure in perovskite thin films. 

Turning now to our results for the PbTiOa (011) sur- 
faces, we find that all metal atoms in the outermost 
layer move inwards irrespective of the termination. Sur- 
face oxygen atoms are displaced outwards for the TiO- 
terminated surface, while oxygen atoms move inwards in 
the O-terminated surface. The displacement patterns of 
atoms in the outermost surface layers are similar to those 
of the BaTi03 (011) surfaces, as well as classical shell 
model results for BaTiOsJ^ For example, the atomic 
displacement magnitudes of Ti and oxygen atoms in 
the TiO-terminated PbTi0 3 (011) surface are -0.0813 a 
and 0.033 ao respectively. The Pb atom is displaced in- 
wards by 0.1194ao for the Pb-terminated surface, sim- 
ilar to the corresponding BaTi03 case. Overall, Table 
III shows similar displacement patterns for the BaTi03 
and PbTiOs (011) surfaces, as well as qualitatively sim- 
ilar results for both ab initio and classical shell-model 
descriptions. 



C. BaTi0 3 and PbTi0 3 (001) and (Oil) surface 
energies 

In the present work, we define the unrelaxed surface 
energy of a given surface termination X to be one half 
of the energy needed to cleave the crystal rigidly into an 
unrelaxed surface X and an unrelaxed surface with the 
complementary termination X'. For BaTi03, for exam- 
ple, the unrelaxed surface energies of the complementary 
BaO and Ti02-terminated (001) surfaces are equal, as 
are those of the TiO and Ba-terminated (011) surfaces 
(and similarly for PbTiOs). The relaxed surface energy 
is defined to be the energy of the unrelaxed surface plus 
the (negative) surface relaxation energy. These defini- 
tions are chosen for consistency with Refs. [l7ll2"8j . Unlike 
the authors of Refs. [27.29], we have made no effort to 
introduce chemical potentials here, so the results must 
be used with caution when addressing questions of the 
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TABLE III: Calculated surface relaxations of BaTi03 and PbTiOa (Oil) surfaces (in percent of the lattice constant) for the 
three surface terminations. 'SM' indicates comparative results from the shell-model calculation of Ref. [lul ]. 
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relative stability of surfaces with different stoichiometry. 

With these definitions, and using the 7-layer slab ge- 
ometries specified in Sec. II B, the energy of the unrelaxed 
BaTi0 3 (001) surface is 

iC r f (X) = \[EZU^O) + £ slab (Ti0 2 ) - 7£ bulk ] (1) 

where X = BaO or TiC>2 specifies the termination, 
JS^(BaO) and £^ b (Ti0 2 ) are the unrelaxed BaO and 
Ti02-terminated slab energies, -Bbuik is energy per bulk 
BaTiC>3 unit cell, and the factor of four comes from the 
fact that four surfaces are created by the two cleavages 
needed to make the two slabs. The relaxation energy for 
each termination can be computed from the correspond- 
ing slab alone using 

A< rf (X) = l[E slah (X) - EZl(X)}, (2) 

where E B iab{X) is a slab energy after relaxation. The 
relaxed surface energy is then 

E SU1 ,(X) = E™? f (X) + ae::Ux). (3) 

Similarly, for the BaTi0 3 (011) cleavage on a 

bulk BaTiO plane gives rise to the complementary TiO 
and Ba-terminated surfaces shown in Fig. 2(c) and (d) 
respectively. Thus, 

EZ r i (X) = ~[£Xb(Ba) + £Xb(TiO) - 6£ bulk ], (4) 

where the energy is the same for X = TiO or Ba, 
£ SSb 1 ( Ba ) and £ s u iab el ( T iO) are energies of the unrelaxed 
slabs. Relaxation energies can again be computed inde- 
pendently for each slab in a manner similar to Eq. (2). 



Finally, the (011) surface can also be cleaved to give 
two identical self-complementary O-terminated surfaces 
of the kind shown in Fig. 2(b). In this case the 7-layer 
slab has the stoichiometry of three bulk unit cells, so the 
relaxed surface energy of the O-terminated (011) surface 
is 

E suii (0) = i[£ slab (0) - 3£ bulk ], (5) 

where E s \ a ^(0) is the relaxed energy of the slab having 
two O-terminated surfaces. Everything said here about 
BaTiOa surfaces applies in exactly the same way to the 
corresponding PbTi03 surfaces. 

The calculated surface energies of the relaxed BaTiOs 
(001) and (011) surfaces are presented in Table IV. In 
BaTi03, the relaxation energies of the Ti02 and BaO- 



TABLE IV: Calculated surface energies for BaTiOg and 
PbTi0 3 (001) and (011) surfaces (in eV per surface cell). 'SM' 
indicates comparative results from the shell-model calculation 
of Ref. [l||. 
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TABLE V: Calculated magnitudes of atomic displacements D (in A), effective atomic charges Q (in e), and bond populations 
P between metal-oxygen nearest neighbors (in 10 -3 e) for the BaTiOa and PbTi03 (001) surfaces. 

BaTiOa (001) surface PbTi0 3 (001) surface 
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D 


O 


—0.014 


O 


—0.025 


O 


0.012 


O 


—0.012 


Q 




— 1.280 




—1.473 




— 1.184 




—1.128 


P 




—38 




80 




44 




106 


D 


Ba 


0.101 


Ti 


0.070 


Pb 


0.209 


Ti 


0.121 


Q 




1.767 




2.379 




1.275 




2.331 


p 




—30 




88 




8 




80 


D 


o 


0.015 


O 


0.056 


O 


0.050 


O 


0.091 


Q 




-1.343 




— 1.418 




-1.167 




-1.258 


p 




90 




-30 




80 




6 


Q 


Ti 


2.365 


Ba 


1.803 


Ti 


2.335 


Pb 


1.358 


p 




104 




-36 




108 




24 


Q 


O 


-1.371 


O 


-1.417 


O 


-1.207 





-1.259 


p 




-34 




98 




18 




96 


Q 


Ba 


1.797 


Ba 


1.797 


Pb 


1.354 


Pb 


1.354 


p 




-34 




-34 




16 




16 


Q 


O 


-1.388 


O 


-1.388 





-1.232 





-1.232 


p 




98 




98 




98 




98 


Q 


Ti 


2.367 


Ti 


2.367 


Ti 


2.341 


Ti 


2.341 



terminated surfaces (—0.23 and — O.lleV respectively) 
are comparable, leading to rather similar surface ener- 
gies. On the (011) surfaces, however, the relaxation 
energies vary more strongly with termination. For ex- 
ample, we find a relaxation energy of — 2.13eV for the 
TiO-tcrminated surface, much larger than — 0.93eV for 
the Ba-terminated surface. The relaxation energy of 
— 1.15eV for O-terminated surface gives rise to a relaxed 
energy of the O-terminated surface (1.72 eV) that is much 
lower than the average of the TiO and Ba-terminated sur- 
faces (2.64 eV), indicating that it takes much less energy 
to cleave on an O2 plane than on a BaTiO plane. The 
shell-model results of Ref. [l6[ for the BaTi03 surfaces 
are given shown for comparison; the results are qualita- 
tively similar, but there are some significant quantitative 
differences, especially for the Ba-terminated (011) sur- 
face. 

The corresponding results are also given for the (001) 
and (011) surfaces of PbTi0 3 in Table IV. The results 
for the (001) surfaces are similar to those for BaTiOs, al- 
though the relaxed surface energies are somewhat lower. 
For the case of the (011) surfaces, however, we find a dif- 
ferent pattern than for BaTi03. We find a very large 
relaxation energy of — 1.75eV for the TiO-terminated 
surface, compared with — 1.08eV for the Pb-terminated 
surface and — 1.12eV for the O-terminated surface. The 
average energy of the TiO and Pb-terminated surfaces 
is now 1.69 eV, to be compared with 1.72 eV for the O- 
terminated surface, indicating that the cleavage on a Pb- 
TiO or an O2 plane has almost exactly the same energy 
cost. 



D. BaTi0 3 and PbTi0 3 (001) and (Oil) surface 
charge distribution and chemical bonding 

To characterize the chemical bonding and covalency 
effects, we used a standard Mulliken population analysis 
for the effective static atomic charges Q and other local 
properties of the electronic structure as described, for ex- 
ample, in Ref. 3!). hj. The results are presented in Table 
V. Our calculated effective charges for bulk PbTiOa are 
+1.354e for the Pb atom, +2.341 e for the Ti atom, and 
— 1.232 e for the O atom. The bond population describing 
the chemical bonding is +98 me between Ti and O atoms, 
+16 me between Pb and O atoms, and +2 me between Pb 
and Ti atoms. Our calculated effective charges for bulk 
BaTi0 3 are +1.797e for the Ba atom, +2.367e for the 
Ti atom, and — 1.388 e for the O atom indicate a high 
degree of BaTiOa chemical bond covalency. The bond 
population between Ti and O atoms in BaTiOs bulk is 
exactly the same as in PbTiOs, while that between Ba 
and Ti is slightly negative, suggesting a repulsive inter- 
action between these atoms in the bulk of the BaTiOs 
crystal. 

For the Ti0 2 -terminated BaTi0 3 and PbTi0 3 (001) 
surfaces, the major effect observed here is a strength- 
ening of the Ti-0 chemical bond near the BaTiOs and 
PbTiOs (001) surfaces, which was already pronounced 
for the both materials in the bulk. Note that the Ti 
and O effective charges for bulk BaTiOs and PbTiOs 
are much smaller than those expected in an ionic model 
(+4e, and — 2e), and that the Ti-0 chemical bonds in 
bulk BaTiOs and PbTiOs are fairly heavily populated 
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TABLE VI: The A-B bond populations P (in lCT 3 e) and 
interatomic distances R (in A) on (Oil) surfaces of BaTiOa 
and PbTi0 3 . Symbols I-IV denote the number of each plane 
enumerated from the surface. The nearest-neighbor Ti-0 dis- 
tance is 2.004 A and 1.968 A in bulk BaTi0 3 and PbTi0 3 , 
respectively. 



BaTiOs (Oil) surface PbTiOs (Oil) surface 



Atom A 


Atom B 


P R Atom A 


Atom B 


p 


R 


TiO-tcrminated 


TiO-terminated 






Ti(I) 


0(1) 


130 2.047 Ti(I) 


0(1) 


132 


2.019 




O(II) 


198 1.784 


O(II) 


196 


1.766 


O(II) 


Ti(III) 


112 2.009 O(II) 


Ti(III) 


120 


1.948 




Ba(III) 


-24 2.808 


Pb(III) 


24 


2.826 




O(III) 


-26 2.837 


O(III) 


— 20 


2.857 


Ti(III) 


Ba(III) 


-2 3.471 Ti(III) 


Pb(III) 


2 


3.410 




O(III) 


118 2.004 


O(III) 


108 


1.975 




O(IV) 


96 2.004 


O(IV) 


88 


1.976 


Ba(III) 


O(III) 


-32 2.834 Pb(III) 


O(III) 


20 


2.783 




O(IV) 


-38 2.816 


O(IV) 


8 


2.734 


O(III) 


O(IV) 


-30 2.834 0(111) 


O(IV) 


-36 


2.706 


Ba-terminated 


Pb-terminated 






Ba(I) 


O(II) 


-38 2.664 Pb(I) 


O(II) 


126 


2.589 


O(II) 


Ba(III) 


-36 2.850 O(II) 


Pb(III) 


24 


2.742 




Ti(III) 


84 2.022 


Ti(III) 


74 


1.902 




O(III) 


-38 2.859 


O(III) 


—46 


2.739 


Ba(III) 


O(III) 


-36 2.834 Pb(III) 


O(III) 


— 10 


2.783 




O(IV) 


-36 2.834 


O(IV) 


36 


2.813 


Ti(III) 


O(III) 


76 2.004 Ti(III) 


O(III) 


62 


1.968 




Ba(III) 


-2 3.471 


Pb(III) 





3.408 




O(IV) 


98 2.008 


O(IV) 


92 


2.018 


O(III) 


O(IV) 


-46 2.825 O(III) 


O(IV) 


-46 


2.816 


O-terminated 


O-terminated 






0(1) 


Ba(II) 


-26 2.697 0(1) 


Pb(II) 


50 


2.503 




Ti(II) 


168 1.722 


Ti(II) 


128 


1.694 




O(II) 


-24 2.801 


O(II) 


-26 


2.689 


Ba(II) 


O(II) 


-40 2.664 Pb(II) 


O(II) 


78 


2.574 




Ti(II) 


-2 3.306 


Ti(II) 


4 


3.094 


Ti(II) 


O(II) 


82 2.040 Ti(II) 


O(II) 


92 


1.977 




O(III) 


112 1.689 


O(III) 


126 


1.831 


O(II) 


O(III) 


-12 2.825 O(II) 


O(III) 


-24 


2.779 


Ba(II) 


O(III) 


-10 2.968 Pb(II) 


O(III) 


26 


2.716 


O(III) 


O(IV) 


-14 2.945 O(III) 


O(IV) 


-42 


2.842 




Ti(IV) 


60 2.159 


Ti(IV) 


60 


2.051 




Ba(IV) 


-24 2.767 


Pb(IV) 


-2 


2.712 



for both materials. 

The Ti-0 bond population for the Ti02-tcrminated 
BaTi0 3 and PbTi0 3 (001) surfaces are +126 me and 
+114 me respectively, which is about 20% larger than the 
relevant value in the bulk. In contrast, the Pb-0 bond 
population of +54 me) is small for the PbO-terminatcd 
PbTiC>3 (001) surface, and the Ba-0 bond population of 
—30 me is even negative for the BaO-terminated BaTiC>3 
(001) surface, indicating a repulsive character. The effect 
of the difference in the chemical bonding is also well seen 
from the Pb and Ba effective charges in the first surface 
layer, which are close to the formal ionic charge of +2 e 
only in the case of the BaTiOs crystal. 

The interatomic bond populations for three possible 



TABLE VII: Calculated Mulliken atomic charges Q (in e) 
and changes in atomic charges AQ with respect to the bulk 
charges (in e) on (011) surfaces of BaTiOs and PbTi0 3 . for 
three terminations. The Mulliken charges are 2.341 e for Ti, 
-1.232 e for O, and 1.354 e for Pb in bulk PbTi0 3 , and 2.367 e 
for Ti, -1.388 e for O, and 1.797 e for Ba in bulk BaTi0 3 . 

BaTiOs (011) surface PbTi0 3 (011) surface 



Atom 


Q 


AQ 


Atom 


Q 


AQ 


TiO-terminated 




TiO-terminated 




Ti(I) 


2.216 


-0.151 


Ti(I) 


2.212 


-0.129 


O(I) 


-1.316 


0.072 


O(I) 


-1.261 


-0.029 


O(II) 


-1.155 


0.233 


O(II) 


-1.057 


0.175 


Ba(III) 


1.757 


-0.04 


Pb(III) 


1.253 


-0.101 


Ti(III) 


2.353 


-0.014 


Ti(III) 


2.328 


-0.013 


O(III) 


-1.299 


0.089 


0(111) 


-1.18 


0.052 


O(IV) 


-1.402 


-0.014 


O(IV) 


-1.239 


-0.007 


Ba— terminated 




Pb— terminated 




Ba(I) 


1.636 


-0.161 


Pb(I) 


1.122 


-0.232 


O(II) 


-1.483 


-0.095 


O(II) 


-1.140 


0.092 


Ba(III) 


1.799 


0.002 


Pb(III) 


1.340 


-0.014 


Ti(III) 


2.368 


0.001 


Ti(III) 


2.343 


0.002 


O(III) 


-1.446 


-0.058 


O(III) 


-1.277 


-0.045 


O(IV) 


-1.392 


-0.004 


O(IV) 


-1.247 


-0.015 


O— terminated 




O—terminated 




0(1) 


-1.158 


0.23 


O(I) 


-1.011 


0.221 


Ba(II) 


1.766 


-0.031 


Pb(II) 


1.257 


-0.097 


Ti(II) 


2.213 


-0.154 


Ti(II) 


2.237 


-0.104 


O(II) 


-1.452 


-0.064 


O(II) 


-1.261 


-0.029 


O(III) 


-1.317 


0.071 


0(111) 


-1.215 


0.017 


Ba(IV) 


1.792 


-0.005 


Pb(IV) 


1.355 


0.001 


Ti(IV) 


2.317 


-0.05 


Ti(IV) 


2.317 


-0.024 


O(IV) 


-1.407 


-0.019 


O(IV) 


-1.233 


-0.001 



BaTiC>3 and PbTiC>3 (011) surface terminations are given 
in Table VI. The major effect observed here is a strong in- 
crease of the Ti-0 chemical bonding near the BaTiC>3 and 
PbTiC>3 (Oil) surface as compared to already large bond- 
ing near the (001) surface (+126 me and +114me, respec- 
tively) and in the bulk (+98 me). For the O-terminated 
(Oil) surface the 0(I)-Ti(II) bond population is as large 
as +168 me for BaTi0 3 and +128 me for PbTiOg, i.e., 
considerably larger than in the bulk and on the (001) 
surface. 

Our calculations demonstrate that for the TiO- 
terminated BaTi0 3 and PbTi0 3 (011) surfaces, the Ti-0 
bond populations are larger in the direction perpendic- 
ular to the surface (+198 me for BaTi0 3 and +196 me 
for PbTi0 3 ) than in plane (+130 me for BaTiOs and 
+132 me for PbTi0 3 ). The Ti-0 bond populations for 
the TiO-terminated PbTi0 3 (Oil) surface in the direc- 
tion perpendicular to the surface is twice as large as the 
Ti-0 bond population in PbTi0 3 bulk. 

In Table VII we present the calculated Mulliken effec- 
tive charges Q, and their changes AQ with respect to the 
bulk values, near the surface. We analyzed the charge 
redistribution between different layers in slabs with all 
three BaTi0 3 and PbTi03 (011) surface terminations. 
The charge of the surface Ti atoms in the TiO-terminated 
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BaTi0 3 and PbTi0 3 (001) surface is reduced by 0.151 e 
and 0.129 e, respectively. Metal atoms in the third layer 
lose much less charge. Except in the central layer (and, 
in the case of PbTiOs, in the subsurface layer), the O 
ions also reduce their charges, becoming less negative. 
The largest charge change is observed for BaTiC>3 and 
PbTiC>3 subsurface O atoms (+0.233 e and +0.175e, re- 
spectively). This gives a large positive change of +0.466 e 
and +0.350 e in the charge for each BaTiC>3 and PbTiC>3 
subsurface layer. 

On the Ba-terminated and Pb-terminated BaTiC>3 and 
PbTiC>3 (011) surface, negative changes in the charge 
are observed for all atoms except for Ba and Ti in the 
BaTiC>3 third layer, Ti atom in the PbTiC>3 third layer, 
and subsurface oxygen atom in PbTiC>3. The largest 
charge changes are at the surface Ba and Pb ions. It 
is interesting to notice that, due to the tiny difference in 
the chemical bonding between BaTiOs and PbTiOs per- 
ovskites, the charge change for the BaTiOs subsurface O 
ion (-0.095 e) and PbTiOg subsurface O ion (+0.092 e) 
have practically the same magnitude, but opposite signs. 

For the O-terminated BaTi0 3 and PbTi0 3 (011) sur- 
faces, the largest calculated changes in the charge are 
observed for the BaTi0 3 and PbTiOs surface O atom 
(+0.230 e and +0.221 e, respectively). The change of the 
total charge in the second layer is negative and almost 
equal for both materials. For the BaTiOs crystal, this re- 
duction by 0.249 e comes mostly from Ti atom (— 0.154 e). 
In the PbTiC>3 crystal, the reduction by 0.230 e appears 
mostly due to a decrease of the Ti atom charge by 0.104 e, 
as well as a reduction of the Pb atom charge by 0.097 e. 

IV. CONCLUSIONS 

In summary, motivated by the scarcity of experimen- 
tal investigations of the BaTiOs and PbTiOs surfaces 
and the contradictory experimental results obtained for 
the related SrTiC>3 surface ) 19 i 20 we have carried out pre- 
dictive electronic structure calculation to investigate the 
surface atomic and electronic structure of the BaTiOs 
and PbTiOs (001) and (011) surfaces. Using a hybrid 
B3PW approach, we have calculated the surface relax- 



ation of the two possible terminations (TiC>2 and BaO 
or PbO) of the BaTi0 3 and PbTi0 3 (001) surfaces, and 
three possible terminations (TiO, Ba or Pb, and O) of the 
BaTi0 3 and PbTi0 3 (011) surfaces. The data obtained 
for the surface structures are in a good agreement with 
previous LDA calculations of Padilla and Vanderbilt, 4 
the LDA plane- wave calculations of Meyer et al.^ and 
in fair agreement with the shell-model calculations of 
Heifets et al^ 

According to our calculations, atoms of the first sur- 
face layer relax inwards for BaO and PbO terminated 
(001) surfaces of both materials. Outward relaxations 
of all atoms in the second layer are found at both ter- 
minations of BaTiOs and PbTiOs (001) surfaces. In 
BaTi03, the rumpling of the Ti02-terminated (001) sur- 
face is predicted to exceed that of the BaO-terminated 
(001) surface by a factor of two. In contrast, PbTiOs 
exhibits practically the same rumplings for both (Ti02 
and PbO) terminations. Our calculated surface energies 
show that the Ti02-terminated (001) surface is slightly 
more stable for both materials than the BaO or PbO- 
terminated (001) surface. The O-terminated BaTiOs and 
TiO-terminated PbTi03 (011) surfaces have surface en- 
ergies close to that of the (001) surface. Our calcula- 
tions suggest that the most unfavorable (011) surfaces are 
the Ba or Pb-terminated ones for both the BaTiOs and 
PbTi03 cases. We found that relaxation of the BaTiOs 
and PbTiOs surfaces for is considerably stronger for all 
three (011) terminations than for the (001) surfaces. The 
atomic displacements in the third plane from the surface 
for the three terminations of BaTiOs and PbTiOs (011) 
surfaces are still large. Finally, our ab initio calculations 
indicate a considerable increase of Ti-0 bond covalency 
near the BaTiOs and PbTiOs (011) surface relative to 
BaTi0 3 and PbTi0 3 bulk, much larger than for the (001) 
surface. 
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